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Abstract 


A two-state Cellular Automaton (CA) with neighborhood N is 
called an Interval Life-Like Freezing Cellular if the local rule f can 
be defined from an interval Z C [0,...,|N|], such that: (1) when a cell 
c is in state 0, it becomes 1 if and only if the number of neighbors of 
c in state 1 belongs to Z. (2) if c is in 1, then it remains in state 1 
independently on the states of their neighbors. 

Recently (Goles et al. Information and Computation 2020), an 
exhaustive study of the dynamics and complexity of ILLFCA was car- 
ried out in the von Neumann neighborhood, where the complexity and 
dynamics of such CAs were classified in five groups, namely Trivial, 
Algebraic, Topological, Fractal and P-Complete rules. 

In this exploratory paper, we propose a follow-up of the later re- 
search, extending it to ILLFCA equipped with the Moore neighbor- 
hood. Interestingly, this class contains the famous Life-without-Death, 
also known as the freezing version of Conway Game of Life. In fact, we 
show that the family of ILLFCA with Moore neighborhood has a rich 
diversity of rules, some of them not classifiable in the previous groups. 
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1 Introduction 


In this exploratory paper, our goal is to explore an exhaustive classification 
of the complexity of the Interval Life-Like Freezing Cellular Automata with 
Moore neighborhood. We begin defining this class. 

A freezing cellular automaton is a CA that admits an order on its states, 
such that the state evolution of any node is non-decreasing in any orbit. 
Several models that received a lot of attention in the literature are ac- 
tually freezing dynamics, for instance: self-assembly tilings [1], forest fire 
or epidemic [2, 3] propagation models, bootstrap percolation [4], etc. On 
the other hand, their complexity as computational models has been stud- 
ied from various points of view: as language recognizers where they corre- 
spond to bounded change or bounded communication models [5], for their 
computational universality [6, 7], as well as for various associated decision 
problems [8]. 

A Life-Like Cellular Automata are a family of CAs that generalize Con- 
way’s Game of Life [9]. This family contains the totalistic CA with state 
{0,1} and Moore neighborhood N such that there exist sets B, S C {0,1..., 8} 
such that: 


flen) =1s (20,0) =0A >» Li € B) V (x(0,0) =LA `> Li € S) 
zEN zEN 


Born rule Survivor rule 


A Life-Like CA is called Interval Life-Like Cellular Automaton (ILLFCA) if 
both sets B and S are closed intervals. For instance, the Game of Life is the 
Interval Life-Like CA with B = {3}, S = {2,3}. Observe that an Interval 
Life-Like CA is freezing if S = {0,1...,8}. The freezing version of the Game 
of Life is known as the Life Without Death [10] (B = {3}, S = {0,1,...,8}). 

Since we only consider Interval Life-Like CA that are freezing, we omit 
the set S in their definition, as these rules are only defined by the interval 
B. Moreover, for 0 < a < b < 8, we denote Ra» the ILLFCA such that 
B = [a,b]. For instance the Life Without Death is denoted R33. In order 
to illustrate the different rules that belong to this family, we present them 
in the table shown in Figure 1. In the i-th row of the latter table we put all 
the ILLFCAs whose interval starts at à and in column j all those ILLFCAs 
whose interval closes at 7. Thus, we obtain an upper triangular matrix (see 
Figure 1). From here we identify four types of rules: 


e The rules Rig for i € {1,...,7}, called threshold rules. In these rules 
a cell switches to state 1 if at least a threshold à of its neighbors are 
in state 1. These rules are denoted by 6;. 


e The rules Ry; for j € {2,...,7}, called anti-threshold rules. In these 
rules a cell switches to state 1 if the set of neighbors in state 1 is at 
most i. These rules are denoted 6;. 


e Finally, the rules R;; for i € {1,...8}, where a cell is born if there 
are exactly à living neighbors are called delta rules. These rules are 
denoted 6;. 


e The inner rules are all the remaining rules, that is to say, all rules Rij 
such that 1<i< j <8. 
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Figure 1: Family of ILLFCAs. In blue we highlight the threshold rules, in 


gray the anti-threshold rules, and in red the delta rules. All the remaining 
are inner rules. 


We remark that we do not study the rules Roy. This is because all these 
rules satisfy the following property: after the first iteration, all cells with 
0 alive neighbors become alive. Then, we can consider that in the next 
iterations of the dynamics the rule R1, is applied. Thus, studying the Ro» 
rules is equivalent to study the R1,4 rules with initial configurations without 
isolated neighbors. 

In order to measure the complexity of each rule, we consider the com- 
putational complexity of the decision problem consisting in predicting the 
future state of each cell. More precisely, given an initial condition and a cell 
in state 0 we would like to predict if it will eventually change to state 1 or it 
will be stay fixed in state 0. Observe that since the rules we are considering 
are freezing, this latter task corresponds to study of the fixed points [11]. We 
say that a cell is stable, if it is initially not alive and never changes. Thus, 
when the initial configuration reaches a fixed point, the remaining non-living 
cells are the stable cells. Conversely, the alive cells at the fixed point are 
the non-stable (unstable) cells or they are cells that were initially in state 1. 
The definition of stability introduces the following decision problem: 


Stability r, , 

Input: a periodic configuration x described by {0, 1}]x I] 
and a site u € [n] x [n]. 

Question: Is u a stable cell for x? 


The computational complexity of a decision problem can be defined as 
the amount of resources (i.e. time and space) required to solve it. Observe 
that, since the rules are freezing, a fixed point is always reached in poly- 
nomial time in the size of the initial condition. Indeed, in every time-step 
before reaching an attractor at least one cell is frozen in state 1. Thus, there 
is always a (deterministic) polynomial-time algorithm for solving Stability, 
which consists in simulating the rule Rap- Considering this latter obser- 
vation (and from a classical computational complexity standpoint) we can 
roughly identify the following well-known time-complexity classes: P, the 
class of problems solvable in polynomial time on a sequential computer, and 
NC, the class of problems solvable in poly-logarithmic time in a PRAM ma- 
chine, with a polynomial number of processors [12]. Roughly speaking, NC 
is the class of problems for which there exists a fast parallel algorithm which 
solves them. Observe that, since any parallel algorithm can be straightfor- 
wardly executed sequentially (the execution times depends on the number of 
processors and the time that each processor takes in the parallel computa- 
tion), we have that NC C P. In this context, it is well-believed that NC # 
P, and so, there exist “inherently sequential” problems, i.e. problems that 
are in P but they are not in NC. Those inherently sequential problems are 
usually called P-complete problems. More precisely, a decision problem is 
P-complete if it is in P and if any other problem in P is reducible (by an 
NC-reduction or a logarithmic space reduction) to that decision problem. 
Thus, if for some P-complete problem, there is a fast parallel algorithm 
solving it, then necessarily P=NC [12, 13]. 

Griffeath and Moore [14] showed that Stability is P-complete for the Life 
Without Death. This was accomplished by simulating arbitrary boolean 
circuits by a series of gadgets representing logic gates as configuration pat- 
terns. This was similar to the work done by Banks [15]. On the other hand, 
Goles et al. [8] showed that on a regular graph of degree n, the stability 
problem for a threshold rule with threshold n — 1 is in NC. Thus, for 67, 
Stability is in NC. 

Considering the von Neumann neighborhood Goles et al. [16] studied 
Stability for two-state freezing totalistic automata. They classified the latter 
rules in five groups: 


Trivial: It is easy to determine whether a cell will change. Dynamics are 
of the following types: any initial pattern fills everything, initial con- 
figurations remain unchanged, or a fixed point is reached in a few 
iterations. 


Algebraic: The rule, or part of it, can be computed as a composition of 
associative operation, which is used to quickly compute the fixed point 
of any initial configuration. Dynamics are waves of live cells that 
collide and, depending on the position of the cells at the point of 
collision, may or may not change. 


Topological: The set of stable cells can be characterized according to the 
connectedness of the induced graph in the grid. In the case of the 
strict majority, the characterization is that the stable cells are initially 
inactive cells that form biconnected components or paths that connect 
them. In the case of the non-strict majority, the stable cells are those 
that form a 3-connected component of initially inactive cells. 


P-complete: are rules capable of simulating logic circuits. This is the only 
case where Stability was shown to be P-complete. 


Fractal: The dynamics are of patterns that grow in a self-repeating way. 
For this set, the complexity of Stability was not settled. 


2 A landscape of Interval Life-like Freezing CA 


By analyzing the rules and performing numerical experiments, we observe 
that the ILLFCA with the Moore neighborhood can be grouped in a similar 
way than it was done in [16] for the von Neumann case. Nevertheless, we 
find a new set of rules, that we denote signal rules, which exhibit a behavior 
not found in the von Neumann case. We now detail the results regarding 
the threshold, anti-threshold and delta rules. 
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Figure 2: Classification according to the observed in computational simula- 
tions. For those marked with *analytical complexity proofs are provided. 


Trivial rules: R13, R23, R33, Ri,7, R77 and Rss. It is easy to see that for 
rule R18 the two possible fixed points are all alive or all dead configurations. 
Rule R; 7 is similar, but admits some fixed points with dead cells having all 
neighbors alive. Rule Rgs is trivial as every configuration reaches a fixed 
point in one time-step. For rule R2 we have that a configuration has at least 
two alive cells sharing the same neighborhood, then everything converges to 
the configurations where all the cells are alive. Finally, the case of the rule 


R3,8, is slightly more complicated as there are three-cell configurations that 
are stable (see Figure 4), but it is enough to study all patterns of a given 
fixed size. To the following figures, white means death and black means 
alive. 


Figure 3: Pattern that forces any Figure 4: Stable pattern for rule 

configuration to reach the all alive R338. 

fixed point for rule R338. 
Fractal rules: R11, R12, R22. In Figures 5, 6 we show examples of fractal 
dynamical behavior exhibited by rules R, and R12 from a single alive 
cell. Rule R22 requires more than one cell to exhibit such a behavior, since 
otherwise it does not change. In Figure 7 it is shown the behavior starting 
from three initial alive cells. 


Figure 5: Rule R11 Figure 6: Rule R1 2 Figure 7: Rule Ro 2 
Algebraic rules: R13, R14, R15, R16. In Figure 8 we show an example of 
the algebraic dynamic of all these rules. For all these rules, waves starting 
from alive cells propagate. Stable dead cells appear in the interface in which 
the different waves collide. 


Initial Rule R13. Rule HE Rule R15. 
configuration. 

Figure 8: Example of dynamics for the rules R13, R14, and R15 starting 

from the same initial configuration. 


Topological rules: rules R48, Rs5,3, Res, R73, R77. We analyze each of 
these rules separately. 


Rule R48: This is the majority rule (also known as Bootstrap percola- 
tion). It can be shown that stable patterns are necessarily 4-connected. 
In addition, a sufficient condition for stability is 5-connectedness. Nev- 
ertheless, there are stable structures that are not 5-connected. 


Figure 10: An example of 


a stable configuration which 
is not 5-connected. 


Figure 9: An example of a 5- 
connected pattern of stable 
cells. 


Rule Rsg: This is the freezing strict majority rule. Similar to the case 
of the von Neumann neighborhood, there exist stable bounded 4- 
connected components. Non-4-connected structures can be stabilized 
by attaching them to other 4-connected components (see Figures 11 
and 12). 


Figure 11: Stable bounded Figure 12: Two bounded 
configuration. configurations stabilizing a 
third non-4-connected one. 


Rule Re: In this case, there are bounded stable configurations. A suffi- 
cient condition for stability is 3-connectedness. However, there are a 
bounded non 3-connected components that are stable (see Figure 13). 


Figure 13: Left: Stable bounded configuration. Middle: Three Stable 
bounded configuration stabilizing a cell. Right Eight Stable patterns stabi- 
lizing a set of cells. 


Rule R73: As it is shown in [8], a cell is stable if and only if it is in a 
bi-connected component or in a path between them. 


Rule f7,7: This is almost a trivial rule. In order to decide if a cell is stable, 
it suffices to show that the cell belongs to a straight line of dead cells 
surrounded by alive cells (see Figure 14). 


Figure 14: Example of a stable cell for rule R7,7. Two signals 
arrive simultaneously to red cells, then the red pass from six 
to eight neighbors. 


P-complete rules: R33, R44, R5,5 and R66. Rule R33 is the Life Without 
Death, which is shown to be P-Complete in [14]. For the remaining rules, 
to show P-Completeness, we reduce the monotone circuit value problem 
[17] to Stability. To this end, we simulate an arbitrary monotone circuit 
by exhibiting gadgets which simulate the AND, OR, and XOR gates (see 
Figure 15) . The XOR gate will be used to simulate planar crossing (see 
e.g. [16]). 
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Figure 15: AND, OR and XOR gates (left, center, right) for rules R44, R55, 
and Re, (top, middle, bottom), respectively 


Signal rules: R23. We see that the rule is capable of producing different 
types of signals (see Figures 16 and 17). This behavior could have two 
different implications: on the one hand, it could imply that the rule is 
capable of universal computation, as the case of rule 110 [18]. On the other 
hand, if the interactions of the signals are quite limited, it is possible that 
one can take advantage of that simplicity to develop an efficient algorithm 
for solving Stability. Interestingly, this is a type of behavior which was not 
observed in the von Neumann case. 


Figure 16: Rule R23. Initial Figure 17: Rule R23. Final 
configuration. configuration. 
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Appendix 


Definition 2.1. k-connectivity A connected graph is k-connected if and 
only if removing any k — 1 vertices, the graph still connected. 
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